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The cosmological reconstruction scheme for modified F{R) gravity is developed in terms of e- 
folding (or, redshift). It is demonstrated how any FRW cosmology may emerge from specific F{R) 
theory. The specific examples of well-known cosmological evolution are reconstructed, including 
ACDM cosmology, deceleration with transition to phantom superacceleration era which may develop 
singularity or be transient. The application of this scheme to viable F{R) gravities unifying infiation 
• with dark energy era is proposed. The additional reconstruction of such models leads to non- 

I leading gravitational correction mainly relevant at the early /late universe and helping to pass the 

, cosmological bounds (if necessary). It is also shown how cosmological reconstruction scheme may 

' be generalized in the presence of scalar field. 
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INTRODUCTION 



Modified gravity approach suggests the gravitational alternative for unified description of inflation, dark energy 
and dark matter without the need to introduce by hands the inflaton and extra dark components. Moreover, the easy 
explanation of inflationary or dark energy phase in such scenario follows: the corresponding era is emerging due to 
dominance of the speciflc gravitational sector in the course of the universe expansion. In other words, the early-time 
D '. and late-time acceleration is governed by the universe expansion within the specific modified gravity theory. Special 
interest in this gravitational paradigm for the description of the universe evolution is related with F{R) gravity (for a 
■ general review, see [lj|) due to its quite simple structure if compare with more general modified gravity which includes 
I . all curvature invariants as well as non-local terms. Nevertheless, even in frames of F{R) gravity the background 
J> ' evolution (due to high non-linearity of the problem) is often non-explicit and/or non-analytic process. Prom another 
On , side, any realistic modified gravity should pass not only the local tests but also the observational cosmological bounds. 
' To comply with cosmological bounds, the reconstruction program in any modified gravity has been developed 
" The cosmological reconstruction of F{R) gravity has been considered in refs.jl, 0, 0, S|- turns out that in most 
. ' cases this reconstruction is done in the presence of the auxiliary scalar which may be excluded at the final step 
50 \ so that any FRW cosmology may be realized within specific reconstructed F{R) gravity. However, the weak point 
of so developed reconstruction scheme is that the final function F{R) represents usually some polynomial in the 
positive/negative powers of scalar curvature. On the same time, the viable models have strongly non-linear structure. 

In the present paper we develop the new scheme for cosmological reconstruction of F{R) gravity in terms of e-folding 
(or, redshift z) so that there is no need to use more complicated formulation with auxiliary scalar @, H, Ql- Using such 
technique the number of examples are presented where F{R) gravity is reconstructed so that it gives the well-known 
5^ ' cosmological evolution: ACDM epoch, deceleration/acceleration epoch which is equivalent to presence of phantom 
\ and non-phantom matter, late-time acceleration with the crossing of phantom-divide line, transient phantom epoch 
and oscillating universe. It is shown that some generalization of such technique for viable F{R) gravity is possible, so 
that local tests are usually satisfied. In this way, modified gravity unifying infiation, radiation/matter dominance and 
dark energy eras may be further reconstructed in the early or in the late universe so that the future evolution may 
be different. This opens the way to non-linear reconstruction of reaHstic F{R) gravity. Moreover, it is demonstrated 
that cosmological reconstruction of viable modified gravity may help in the formulation of non-singular models in 
finite-time future. The reconstruction suggests the way to change some cosmological predictions of the theory in the 
past or in the future so that it becomes easier to pass the available observational data. Finally, we show that our 
method works also for modified gravity with scalar theory and any requested cosmology may be realized within such 
theory too. 
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II. COSMOLOGICAL RECONSTRUCTION OF MODIFIED F{R) GRAVITY 

Let us demonstrate that any FRW cosmology may be realized in specific F{R) gravity. The starting action of the 
F{R) gravity (for general review, see is given by 

J d^X^g (^^^ + /:„,attcr) . (1) 

The field equation corresponding to the first FRW equation is: 

^ + 3 + ij) F'{R) - 18i(4H^H + HHj F"{R) + K^p . (2) 

with R = 6H + 12H^. We now rewrite Eq.(l2]) by using a new variable (which is often called e-folding) instead of the 
cosmological time t, N — In-^. The variable N is related with the redshift z by e~^ = ^ = 1 + z. Since ^ = F[-^ 

and therefore ^ = H^m^ + ^mm^ '^^^ rewrite ^ by 

= + 3{H^+ HH') F'{R) - 18{Uh^H' + {H'f + H^H") F"{R) + n^p . (3) 



Here H' = dH/dN and H" = d'^H/dN^. If the matter energy density p is given by a sum of the fiuid densities with 
constant EoS parameter Wi, we find 

- „-3(l+i",) _ „ + -3il+Wi)N /.^ 

i i 

Let the Hubble rate is given in terms of N via the function g{N) as 

H^g{N)^g{-lnil + z)) . (5) 

Then scalar curvature takes the form: R — 6g'{N)g{N) + 12g{N)'^, which could be solved with respect to TV as 
N = N{R). Then by using (01) and |[5]), one can rewrite (HJ as 



= -18 



<fF\R) 
~dW~ 

+3 (, {N iR)f + g' [N (R)) 5 (i?))) ^ - ^ + E P.oao ^'^^'"'^e-^(i+-)^(«) , (6) 

i 

which constitutes a differential equation for F{R), where the variable is scalar curvature R. Instead of g, if we use 
G{N) =g{NY = the expression ^ could be a little bit simplified: 

<PF{R) f^^,^,,^,, 3^, , .A 



(Ag (N iR)f g' [N (i?)) + g {N {R)f g' (N (R))' + g (TV (R))' g" {N (i?)) 



= -C3G{N{R)){AG'{N{R)) + G''{N{R)))^^+[?,G{N{R)) + \g'{N{R))^ 



dR 



F{R) 



i 

Note that the scalar curvature is given by i? = 'iG' {N) + 12G{N) . Hence, when we find F{R) satisfying the differential 
equation |[6l) or ([7]), such F{R) theory admits the solution ([5]). Hence, such F{R) gravity realizes above cosmological 
solution. 

As an example, we reconstruct the F{R) gravity which reproduces the ACDM-era but without real matter. In the 
Einstein gravity, the FRW equation for the ACDM cosmology is given by 



K 



H' - —H^ + poa-^ = —Hi + poa„h-^^ . (8) 



Here Hq and po are constants. The first term in the r.h.s. corresponds to the cosmological constant and the second 
term to the cold dark matter (CDM). The (effective) cosmological constant A in the present universe is given by 
A = 12Hl. Then one gets 

GiN) = Hi + ^poao , (9) 
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and R = 3G"(iV) + 12G{N) = 12H^ + n^paa^ -^e'^^ , which can be solved with respect to A'' as follows, 



Eq.|[7]) takes the following form: 



^ = -oln ^ , . (10) 



= 3 (i? - 9Fo^) {R - UH^,) ^ - f ii? - 9H^ - . (11) 



(PR \2 V dR 2 

By changing the variable from iZ to a; by a; = — 3, Eq.lfTTjl reduces to the hypergeometric differential equation: 

cP F dF 
= x(l - x)— + {j~{a + l3 + l)x) — - al3F . (12) 
aa;^ ax 

Here 

Solution of (fT2l) is given by Gauss' hypergeometric function F{a, P,^; x): 

F{x) ^ AF{a, P, r, x) + Bx^'^F{a - 7 + 1, /? - 7 + 1, 2 - 7; a;) . (14) 

Here A and B are constant. Thus, we demonstrated that modified F{R) gravity may describe the ACDM epoch 
without the need to introduce the effective cosmological constant. 

As an another example, we reconstruct F{R) gravity reproducing the system with non-phantom matter and phantom 
matter in the Einstein gravity, whose FRW equation is given by 

3 

— H'^^Pqa-^ + ppa". (15) 

Here pq, pp, and c are positive constants. When a is small as in the early universe, the first term in the r.h.s. 
dominates and it behaves as the universe described by the Einstein gravity with a matter whose EoS parameter is 
w — —1 + c/3 > — 1, that is, non-phantom like. On the other hand, when a is large as in the late universe, the second 
term dominates and behaves as a phantom-like matter with w = —1 — c/3 < —1. Then since G{N) = g (iV)^ = H^, 
we find 

G = G,e-=^ + Gpe^^ , Gq='^Pqa^\ Gp = yp^ag . (16) 
Then since R = 3G'(7V) + 12G(A^), 

J, R±^R^-AilU-9c^) 

2 (12 + 3c) ' ^ ' 



when c ^ 4 and 



when c = 4. In the following, just for simplicity, we consider c — 4 case. In the case, the non-phantom matter 
corresponding to the first term in the r.h.s. of l|15p could be radiation with w — 1/3. Then Eq.((7]) in this case is given 
by 

^^ 24GpG, , R \ ^ d^FjR) , 9 f 24GpGq , R \ dF{R) F{R) 
° ^ [~^r- + 24 ) ^^R^ + 2 (v R~ + 24 ) -dR 2" ' ^^^^ 

By changing variable i? to a; by i?^ = —576GpGqX, we can rewrite Eq. lfT9|) as 

, d'^F /3 x\ dF F 
0-^l-^)^^n4 + 4h^-2 ' 
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whose solutions are again given by Gauss' hypergeometric function l|14p with 



7=^, a + P+l = -\, aP=\- (21) 

Let us now study a model where the dominant component is phantom-Hke one. Such kind of system can be easily 
expressed in the standard General Relativity when a phantom fluid is considered, where the FRW equation reads 
H'^i'^) ^ 'T'^ph; Here the subscript ph denotes the phantom nature of the fluid. As the EoS for the fluid is given by 
Pph = WphPph with Wph < —1, by using the conservation equation /9ph + 3-ff(H- Wph)/3ph = 0, the solution for the FRW 
equation H^{t) = ^P-ph is well known, and it yields a{t) = ao(ts — t)~^° , where qq is a constant, Hq = — g^j^p^j^j 
and ts is the so-called Rip time. Then, the solution describes the Universe that ends at the Big Rip singularity in the 
time ts. The same behavior can be achieved in F(R) theory with no need to introduce a phantom fluid. The equation 
([7]) can be solved and the expression for the F{R) that reproduces the solution is reconstructed. The expression for 
the Hubble parameter as a function of the number of e-folds is given by H'^{N) = HqB^^^^". Then, the equation |[7]), 
with no matter contribution, takes the form: 

R^^^+AR^ + BFiR)^0, (22) 

where A = —Ho{l + Hq) and B — Il±|^M This equation is the well known Euler equation whose solution yields 



FiR) = CiR"'+ + C2R"'- , where m± = - — "^"^^^^ '^^ . (23) 

Thus, the phantom dark energy cosmology a{t) — ao{ts — can be also obtained in the frame of F{R) theory and 
no phantom fluid is needed. 

We can consider now the model where the transition to the phantom epoch occurs. It has been pointed out that 
F{R) could behave as an effective cosmological constant, such that its current observed value is well reproduced. One 
can reconstruct the model where late-time acceleration is reproduced by an effective cosmological constant and then 
the phantom barrier is crossed (see ref.[5] for such reconstruction in the presence of auxihary scalar). Such transition, 
which may take place at current time, could be achieved in F{R) gravity. The solution considered can be expressed 
as: 

H^ = hJ-^) + iJo = Hic'"^ + Fo , (24) 

where Hi.Hq and a are positive constants. This solution can be constructed in GR when a cosmological constant and 
a phantom fluid are included. In the present case, the solution can be achieved just by an F{R) function, such 
that the transition from non-phantom to phantom epoch is reproduced. Scalar curvature can be written in terms of 
the number of e-folds again. Then, the equation ^ takes the form: 



x{\ - x)F"{x) + 



6 + m \ 1 



6m / 3m 



Tfl A- 4 

F'{x) —Fix) = , (25) 



where x = (12-ffo ~ R)- The equation l(25|) reduces to the hypergeometric differential equation l|14p . so 

the solution is given, as in some of the examples studied above, by the Gauss' hypergeometric function l|15p . whose 
parameters for this case are given by 

1 ^ 3m + 2 ^m + 4 

7 = -—, a + (3 = , a/3=— , 26 

6171 Zm 2m 

and the obtained F{R) gravity produces the FRW cosmology with the late-time crossing of the phantom barrier in 
the universe evolution. 

Another example with transient phantom behavior in F{R) gravity can be achieved by following the same recon- 
struction described above. In this case, we consider the following Hubble parameter: 



H'{N)=Ho\n{ — ]+Hi=HoN + Hi, (27) 



5 



where Hq and Hi are positive constants. For this model, we have a contribution of an eflective cosmological constant, 
and a term that will produce a superaccelerating phase although no future singularity will take place (compare with 
earlier model Q with transient phantom era) . The solution for the model l(27|) can be expressed as a function of time 

H{t) = ^it~to) (28) 

Then, the Universe is superaccelerating, but as it can be seen from l(28|) . in spite of its phantom nature, no future 
singularity occurs. The differential reconstruction equation can be obtained as 

d'^F(x) , , ,dF{x) , , , 

dx dx 

where we have performed a variable change x = H^N + ffi, and the constant parameters are a-i = Hq, a\ = —Hq 

Hi 



hi = — 2^ and 6o = 2_ffo- The equation (|29|) is a kind of the degenerate hypergeometric equation, whose solutions are 
given by the Kummers series K{a, 6; x): 

nR) = K(-2,-^;^^] . (30) 



2Ho 12 

Hence, such F{R) gravity has cosmological solution with the transient phantom behavior which does not evolve to 
future singularity. 

Let us now consider the case where a future contracting Universe is reconstructed in this kind of models. We study 
a model where the universe is currently accelerating, then the future contraction of the Universe occurs. The following 
solution for the Hubble parameter is considered. 

Hit) = 2Hi{to-t) , (31) 

where Hi and to are positive constants. For this example, the Hubble parameter l(3T|) turns negative for t > tp, 
when the Universe starts to contract itself, while for i <C to, the cosmology is typically ACDM one. Using notations 
Hq = 4i7it§ and Hi — AHi and repeating the above calculation, one gets: 

( ~ Hi 12Hq -3Hi- r\ 
FiR)^K^-SHi,-^;^L—^j . (32) 

Hence, the oscillating cosmology ^3^ that describes the asymptotically contracting Universe with a current accelerated 
epoch can be found in specific F{R) gravity. 

Thus, we explicitly demonstrated that F{R) gravity reconstruction is possible for any cosmology under consideration 
without the need to introduce the auxiliary scalar. However, the obtained modified gravity has typically polynomial 
structure with terms which contain positive and negative powers of curvature as in the first such model unifying the 
early-time inflation and late-time acceleration Q. As a rule such models do not pass all the local gravitational tests. 
Some generalization of above cosmological reconstruction is necessary. 



III. COSMOLOGICAL RECONSTRUCTION OF VIABLE F{R) GRAVITY 

In this section, we show how the cosmological reconstruction may be applied to viable modified gravity which passes 
the local gravitational tests. In this way, the non-linear structure of modified F{R) gravity may be accounted for, 
unlike the previous section where only polynomial F{R) structures may be reconstructed. Let us write F{R) |(T]) in 
the following form: F{R) = Fq(R) + Fi{R). Here we choose Fo{R) as a known function like that of GR or one of 
viable F{R) models introduced in @], or viable F{R) theories unifying infiation with dark energy 0, for example 




/o + /i{(i?-i?o)'"+'+<+^}. 
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Using the procedure similar to the one of second section, one gets the reconstruction equation corresponding to ([7]) 

= -9GiNiR))iiG'iN{R)) + G''iN{R)))^^^^+(^3GiNm + lG'iNiR})^ 

2 



-9G (N (R)) (4G' (N (i?)) + G" (N (R))) + (^3G (N (i?)) + ^G' [N (i?))) 



dR 



i 

The above equation can be regarded as a differential equation for Fi{R). For a given G{N) or g{N) ([5]), if one can 
solve (HI) as F{R) = F{R), we also find the solution of §^ as Fi{R) = F{R) - Fo{R). For example, for G{N) 
by using (fT4|) . we find 

Fi (i?) = AF{a, 7; x) + Bx^-^F{a - 7 + 1, /5 - 7 + 1, 2 - 7; x) - Fq (i?) . (35) 
Here a, /3, 7, and a; are given by x = — 3 and (fT3|) . Using Fo{R) (031) one has 



Fi(i?) = AF(a,/3,7;x) + i?xi-^F(a-7+l,/3-7+l,2-7;x)--^ i? ,„^, ° , . 



(36) 



which describes the asymptotically de Sitter universe. Instead of d (fT3|) . if we choose a, 7, and a; 



3H^ 

as i?^ = —576GpGqX and in (|2T1) . Fi{R) ((36l) shows the asymptotically phantom universe behavior, where H diverges 
in future. 

One may start from Fq{R) given by hypergeometric function (fT4|) with a; = 3go(L+4) (^^^o ~ ^^^^ l(26|) . In 

such a model, there occurs Big Rip singularity. Let F{R) be F{R) again given by hypergeometric function (fT4| with 
a; = - 3 and ([13]): 

F{R) = iF(5, (3, 7; £) + Bi^-'iFia - 7 + 1, /3 - 7 + 1, 2 - 7; i) , 

i = ^-3, 7 = -i,« + ^ = -i, a/3 = -i. (37) 

If we choose F{R) = F(i?), the ACDM model emerges. Then choosing Fi(i?) = F{R)-Fo{R), the Big Rip singularity, 
which occurs in Fo{R) model, does not appear and the universe becomes asymptotically de Sitter space. Hence, the 
reconstruction method suggests the way to create the non-singular modified gravity models 0, 0, [ll|. Of course, it 
should be checked that reconstruction term is not large (or it affects only the very early-time/late-time universe) so 
that the theory passes the local tests as it was before the adding of correction term. 
Gauss' hypergeometric function F(a,/3, 7;x) is defined by 

,38, 

r(a)r(/3) r(7 + n)nl 

Since 



-3m - 2 ± - 20m + 4 „ _ ~ -1±5 

ao,/3o = <0, a,/3=— — , 39 

4m 12 

when R is large, Fi(i?) behaves as Fi{R) ~ ^(3m+2+Vm--20m+4)/4m^ ^p.^^ ^^^^^ ^-^^^^ expression, since the total 
F{R) — Fq{R) + Fi{R) is given by F{R) ijST]) . the Big Rip type singularity does not occur. The asymptotic behavior 
of Fi{R) cancels the large R behavior in Fo{R) suggesting the way to present the non-singular cosmological evolution. 
We now consider the case that H and therefore G oscillate as 

G{N) = Go + Gi sin , (40) 
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with positive constants Go, Gi, and Nq. Let the ampUtude of the oscillation is small but the frequency is large: 

Go » ^ , A^o » 1 . (41) 

When Gi = 0, we obtain de Sitter space, where the scalar curvature is a constant R = 12Go. Writing G{N) as 

G - f - Go , (42) 
o 

by using ([7]), one arrives at general relativity: 

F{R) = Co (i? - 6Go) . (43) 

Instead of (|42| . using an arbitrary function F, if we write 

G = Go + F{R) ~ F(12Go) , (44) 

we obtain a general F{R) gravity, which admits de Sitter space solution. When Gi 7^ 0, under the assumption (|4T1) . 
one may identify F{R) in (|43|) with Fq{R). We now write G{N) and the scalar curvature R as 

G(iV) = f - Go + ^g(iV) , R^l2Go + ^r{N) , (45) 

D JVq iVo 

with adequate functions g{N) and r{N). Then since i? = 6g'{N)g{N) + 12g{NY and from l(4T|) . we find 

f N I N \ N N 

g{N) = (n, sin _ + - cos - j , r{N) = 4^o sin ^ + cos - (46) 



F(i?)-co(i?-6Go + ^/(i?)) , (47) 
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i^i(i?) = ^/(i?) , (48) 



^ = gS-sJ^\+o(%N,\ , (49) 



By assuming 



and identifying 



from ([34l) . one obtains 

n = Gn , , , „ , 

dr \No J \Nq 

which can be solved as 

f{R) - (cos-i r T rv/l - ^2) . (50) 

Then at least perturbatively, one can construct a model which exhibits the oscillation of H. 

Before going further, let us find F{R) equivalent to the Einstein gravity with a perfect fiuid with a constant EoS 
parameter w, where H behaves as 

AH2^p„e-3(-+i) . (51) 

Then 

G(N) = ^e-3(-+i) , R(N) - (1 - 3u;) A^^poe-^l-'+D , (52) 



which could be solved as 



1 R 

N = , In -^r— . (53) 

3{w + l) (l-3w)K2po ^ ^ 
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Therefore Eq.|[7l) has the following form: 

3{1 + w) ^ cPF{R) l + 3w dF{R) F{R) 

l-3w dR^ 2(1 -3w) dR 2 ' ^ ' 

whose solutions are given by a sum of powers of R 

F{R) = + . (55) 

Here F± are constants of integration and n± are given by 



. , 7 + 9w , // 7 + 9w Y 2(1 -3w) 1 



If w > —1/3, the universe is decelerating but if — 1 < w < —1/3, the universe is accelerating as in the quintessence 
scenario. 

By using the solution (fT4|) . which mimics ACDM model, and the solution (|55l) . one may consider the following 
model: 

F{x) = {AF{a, (3, 7; ^) + Bx'"^F{a - 7 + 1, - 7 + 1, 2 - 7; a;)} . ."x^.j^^ + ^+^"^ + ' 

(57) 

Here Ri is a constant which is sufficiently small compared with the curvature Rq in the present universe. On the 
other hand, we choose a positive constant A to be large enough. We also choose F± to be small enough so that only 

a(-^-^) 

the first term dominates when R^ Ri. Note that the factor — j-ji — rn — W\ behaves as step function when 

A is large: 

eH^-^) riwheni?>i?i 
eH^-^) + e-(*-^) " ""^^ = i when R < i?, ' ^''^ 

Then in the early universe and in the present universe, only the first term dominates and the ACDM universe could 

a(J^-i) 

be reproduced. In the future universe where i?, <C the factor — f r° \ rn — V decreases very rapidly and the 

o V^^r~ /+e~ / 

second terms in l(57l) dominate. Then if w > —1/3, the universe decelerates again but if — 1 < w < —1/3, the universe 
will be accelerating as in the quintessence scenario. 

Thus, we explicitly demonstrated that the viable F{R) gravity may be reconstructed so that any requested cosmology 
may be realized after the reconstruction. Moreover, one can use the viable F{R) gravity unifying the early-time 
infiation with late-time acceleration (and manifesting the radiation/matter dominance era between accelerations) and 
passing local tests in such a scheme. The (small) correction term Fi{R) can be always constructed so that it slightly 
corrects (if necessary) the cosmological bounds being relevant only at the very early/late universe. This scenario 
opens the way to extremely reahstic description of the universe evolution in F{R) gravity consistent with local tests 
and cosmological bounds. 

IV. RECONSTRUCTION OF MODIFIED GRAVITY WITH EXTRA SCALAR 

We now consider the reconstruction of F(i?)-gravity coupled with a scalar field, whose action is given by 

S = J d^X^ - ia^<^a^0 - V{(t>) + Anattcr) ■ (59) 

Let us redefine the scalar field as = 4'{f)^ 

S = J d*x^(^^-^d^^d^'ip^V{^)+Cn,,uc?j ■ (60) 



Here 

2 



If (j) only depends on the time-coordinate t or e-folding N , wc may choose ip = t or (p = N. 

Then the equations corresponding to the first and second FRW equations have the following form 



= - 



F{R) 



+ 3{H' + HH') F'{R) - 18 (aH^H' + H' {H'f + H^H"^ F"{R) 



F{R) 
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(3i?2 + HH') F'{R) + 6 (l6H^H' - AH' {H'f - H {H'f - AH'H'H" - H^H'"^ F"{R) 



-36 (aH^H' + H' {H'f + H^H"^ {{H'f + HH" + AHH'^ F"'{R) 
which can be rewritten as 

K'ij{ip) {if'f = ^^-2HH'F'{R) + 6 (-AH^H' + 7H' {H'f + H {H'f + AH'H'H" + H^H'" + 3H^H"'^ F"{1 
-36 (aH^H' + H' {H'f + H^H"^ {{H'f + HH" + AHH'^ F'"{R) 

2K'V{ip) = F{R) + {-6H''-AHH')F'{R) 

+6 {28H^H' - H' {H'f - H {H'f - AH'H'H" + H^H'" - 3H^H"^ F"{R) 

+36 {ah'''H' + H^ {H'f + H^H"^ {{H'f + HH" + AHH'^ F"'{R) 

i 

Then if we consider the model given by an adequate function K = K{ip), 

k'<jj{^) = {-2K{ip)K'{ip)F' {6K{ip)K'{ip) + l2K{ipf) 

+6 {-AK{ipfK'{^) + 7K{ipfK'{^f + K{ip)K'{ipf + AK{ipf K' {ip)K" {^) + K{ipfK'"{ip) 
+3K{<^fK"{p)) F" {6K{^)K'{ip) + 12K{ip)') 

-36 {AK{^fK'{^) + K{^fK'{^f + K{^fK"{^)) {K'{vf + K{^)K"{^) + AK{^)K'{^)) 

xF'" {QK{^)K'{^) + l2K{^f) +Y,^H^ + ^0 Pioao '^^+™-)e-3(i+-)-| , 

2K''V{ip) = F {6K{ip)K'{^) + 12K{^f) + {-6K{ipf - AK{ip)K'{ip)) F' {6K{(p)K'{ip) + 12K{y^)') 

+6 {28K{^fK'{ip) - K{^fK'{ipf - K{^)K'{^f - AK{^fK'{^)K"{^) + K{<ffK'"{^) 

-3K{ipfK"{ip)) F" {6K{^)K'{^) + 12K{pf) 

+36 {AK{^fK'{^) + K{^fK'{pf + K{^fK"{^)) {K'{pf + K{^)K"{^) + AK{^)K'{^)) 
xF'" {6K{^)K'{^) + 12K{^)') + ^ «2 _ p.,a-3(i+-.)e-3(i+-.)^ , 

i 

we find a solution is given by 

H{N) = K{N) , ip = N . 
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Hence, it is demonstrated that reconstruction can be extended to the case when modified gravity couples with some 
scalar field. Note that extra scalar may be necessary in the situation when some of cosmological bounds (for instance, 
cosmological perturbations theory which is extremely complicated in F{R) gravity, for a review, see cannot be 
passed within only modified gravity. So far, this is not the case and modified gravity which passes local tests and 
cosmological bounds is available (sl. M flo| . 

V. DISCUSSION 

In summary, we developed general scheme for cosmological reconstruction of modified F{R) gravity in terms of 
e-folding (or redshift) without use of auxiliary scalar in intermediate calculations. Using this method, it is possible 
to construct the specific modified gravity which contains any requested FRW cosmology. The number of F{R) 
gravity examples is found where the following background evolution may be realized:ACDM epoch, deceleration 
with subsequent transition to effective phantom superacceleration leading to Big Rip singularity, deceleration with 
transition to transient phantom phase without future singularity, oscillating universe. It is important that all these 
cosmologies may be realized only by modified gravity without use of any dark components (cosmological constant, 
phantom, quintessence, etc). 

It is shown that our method may be applied to viable F{R) gravities which pass local tests and unify the early- 
time infiation with late-time acceleration. In this case, the additional reconstruction may be made so that correction 
term is not large and it is relevant only in the very early/very late universe. Hence, the purpose of such additional 
reconstruction is only to improve the cosmological predictions if the original theory does not pass correctly the precise 
observational cosmological bounds. For instance, in this way it is possible to formulate the modified gravity without 
finite-time future singularity. It is also demonstrated that the reconstruction scheme may be generalized for the case 
of modified gravity with scalar field. 

The present reconstruction formulation shows that even if specific realistic modified gravity does not pass correctly 
some cosmological bounds (for instance, does not lead to correct cosmological perturbations structure) it may be 
improved with eventually desirable result. Hence, the successful development of such method adds very strong 
argument in favour of unified gravitational alternative for infiation, dark energy and dark matter. 
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